Regularity Properties for a System of Interacting Bessel Processes 



Sebastian Andres and Max-K. von Renesse 
Abstract 

We study the regularity of a diffusion on a simplex with singular drift and reflecting 
boundary condition which describes a finite system of particles on an interval with 
Coulomb interaction and reflection between nearest neighbors. 

... As our main result we establish the Feller property for the process in both cases of 

I repulsion and attraction. In particular the system can be started from any initial 

(<~^ ' state, including multiple point configurations. Moreover we show that the process 

CN ■ is a Euclidean semi-martingale if and only if the interaction is repulsive. Hence, 

, contrary to classical results about reflecting Brownian motion in smooth domains, in 

the attractive regime a construction via a system of Skorokhod SDEs is impossible. 



(N 



o 



Finally, we establish exponential heat kernel gradient estimates in the repulsive regime. 

The main proof for the attractive case is based on potential theory in Sobolev spaces 
with Muckenhoupt weights. 

1 Introduction and Main Results 

The study of Brownian motion or more general diffusions in Euclidean domains with reflect- 
ing boundary condition is a classical subject in stochastic analysis, with strong connection to 
boundary regularity theory for parabolic PDF. Starting from the early works by e.g. Fukushima 
fn\ and Tanaka [2D] , the field has seen perpetual research activity, c.f. HH] (and e.g. [1] for a 
more comprehensive list of references). Typically the reflecting process can be obtained in two 
ways. Either by solving a system of Skorokhod SDE involving the local time at the boundary 



00 . or via Dirichlet form methods, and under suitable smoothness assumptions on the domain and 



the coefficients both approaches are equivalent. 

In this paper we study a very specific singular case, in which the drift coefficients of the operator 
(3 ' may diverge and where the equivalence of the two approaches breaks down but the process 

0^ . exhibits good spatial regularity nevertheless. Our case corresponds to a Dirichlet form obtained 

' from the closure of the quadratic form 

><■ £if,f)= I Nfl^qidx), fec\Ti) 



on Lp'iQ., dq), for a very specific choice of domain i7 C and measure q{dx) (see below). Similar 
variants were studied under the smoothness assumption q G H^'^{Q) and V(logg) G L^{^,dq), 
p > N in [22] and [11] respectively, where a Skorokhod decomposition for the induced process 
still holds. 

The case treated in this paper is obtained by choosing 

n = T.N = {x eR'^ : < < • • • < < 1} 
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and 

N 



q^idx) = ^ TT(x^+^ - x^)^-idx^ dx^-- - dx^, (1.1) 

P i=0 

where by convention x° = 1, x^+^ = 1, = (r(/3)/(r(/3/(iV + l)))^+i)"i is a normalization 
constant and /? > is a free parameter. 

We study the process (X^) generated by (the L^(S7v, '?Af)-closure of) the (pre-)Dirichlet form 

f^(/,/)=/ \Vf\\x)qN{dx), feC^i^N), 
whose generator C extends the operator {L^ , D{L^)) 

with domain 

D{L^) = C^,„ = {/ e C^{^n) I V/ • = on all (n - 1) -dimensional faces of O^n}, 
and 1/ denoting the outward normal field on dT,]\f. 

On the level of formal Ito calculus {L^ , D{L^)) corresponds to an order preserving dynamics 
for the location of N particles in the unit interval which solves the system of coupled Skorokhod 
SDEs 

d^l = (Ti^-^)(^^-l^-^T^)dt + V2dwl + dll-'-dll, i = l,...,N (1.3) 

iV + 1 \xl — x^ x^ — XtJ 

where {wi\ are independent real Brownian motions and {/*} are the collision local times, i.e. 
satisfying 

di\>Q, 11= l\^,^^.^dr,. (1.4) 

JO 

{X^) may thus be considered as a system of coupled two sided real Bessel processes with uni- 
form Bessel dimension 5 = (5/{N + 1). Similar to the standard real Bessel process BES{5) 
with Bessel dimension 5 < 1, the existence of , even with initial condition initial condition 
Xq = X G Sat, is not trivial, nor are its regularity properties. 



Our motivation for studying this process is its relation to the Wasserstein diffusion, c.f. 
In [3] we showed that the normalized empirical measure of the system (|1.3|) converges to the 
Wasserstein diffusion in the high density regime for ^ oo. Hence the regularity properties 
of {X^) may give an indication of the regularity of the Wasserstein diffusion, although in the 
present paper we treat the case when the dimension N is fixed. Here our results read as follows. 

Theorem 1.1. For any /? > 0, the Dirichlet form £^ generates a Feller process {X.) on Sat, 
i.e. the associated transition semigroup on L'^(T,n, qj\[) defines a strongly continuous contraction 
semigroup on the subspace C(Sjv) equipped with the sup-norm topology. 

Moreover, for (3 > (A^-|-l) the associated heat kernel (Pt) is exponentially smoothing on Lipschitz 
functions, i. e. for i > 

Lip(PJ) < exp(-(^^ - l)kN ■ t) Lip(/) (1.5) 
for all f &Lip[Ti]\}), where > is a universal constant depending only on N. 
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In the proof of Theorem 11.11 for the more difficuh case /3 < {N + 1) we use some locahzation 
arguments to exploit the geometric symmetries of the problem. A crucial ingredient for this 
approach is the following ('Markov-') uniqueness result for the operator {L^ ,Cff^^) which is 
interesting in its own right. 

Proposition 1.2. For (3 < 2(N + 1), there is at most one symmetric strongly continuous 
contraction semigroup {Tt)t>o on L'^{Ti]\f,qj\f) whose generator {C,V{C)) extends (L^,C^g^). 

Hence, together with Theorem ll.il we can state the following existence and uniqueness result. 

Corollary 1.3. The formal system of Skorokhod SDEs defines via the associated martingale 
problem a unique diffusion process which can he started everywhere on the (closed) simplex Sjv- 

As for path regularity we obtain the following characterization. 

Theorem 1.4. For any starting point x G ^n, {^^) S Sjv C is a Euclidean semi- 
martingale if and only if (3/{N + 1) > 1. 

In particular we obtain that a Skorokhod decomposition of the process X!^ is impossible if /? is 
small enough. This is in sharp contrast to all aforementioned previous works. Moreover, again 
due to the uniqueness assertion of Proposition 11.21 the following negative result holds true. 

Corollary 1.5. If P/{N + 1) < 1, the system of equations (ll.3p . (ll.4p is ill-posed, i.e. it admits 
no solution in the the sense of ltd calculus. 

Theorems 11.11 and 11.41 generalize the corresponding classical results for the familiy of standard 
real Bessel processes, which are proved in a completely different manner, c.f. [6| I17j. 

Strategy of proof. While the proof of Theorem 11.41 consists of a straightforward application of a 
regularity criterion for Dirichlet processes by Fukushima jl3j , the proof of Theorem 11.11 is more 
involved, and entirely different methods are used in the two respective cases P < N -\- 1 and 
P > N + 1. Both cases are non-trivial from an analytic point of view due to the degenerating 
coefficients and the Neumann boundary condition. 

By comparison the case /? > {N -\- 1) is much easier since q^ is then log-concave. Using a 
recent powerful result by Ambrosio, Savare and Zambotti |2j on the stability of gradient flows 
for the relative entropy functional on Wasserstein the contraction estimate ()1.5p is established 
by smooth approximation and coupling. 

For P < (A^ + 1) the measure q^ is no longer log-concave and to prove the Feller property we 
proceed in four main steps. The first crucial observation is that the reference measure q^ can 
locally be extended to a measure q on the full Euclidean space which lies in the Muckenhoupt 
class A2, allowing for a rich potential theory. In particular the Poincare inequality and doubling 
condition hold which imply via heat kernel estimates the regularity of the induced process on 
R^. The second step is the probabilistic piece of the localization, which corresponds to stopping 
the M'^-valued process at domain boundaries where the measure q is 'tame' and to show that 
the stopped process is again Feller. To this end we employ a version of the Wiener test for 
degenerate elliptic diffusions by Fabes, Jerrison and Kenig [9]. The third step is to use the 
reflection symmetry of the problem which allows to treat the Neumann boundary condition 
indeed via a reflection of the extended M'^-valued process. The fourth step is to establish 
the Markov uniqueness which is crucial in order to justify the identification of the processes 
after localization, and here we shall again depend on the nice potential theory available in the 
Muckenhoupt class. 
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The partial resemblance of our proof for /? < {N + 1) to the classical work by Bass and Hsu [5] 
is no surprise. However, we did not find any similar work in the probability literature where 
potential theory for Muckenhoupt weights was used so extensively as in our case. We think 
that the general strategy outlined here may be useful also in other, less specific situations. 



2 Dirichlet Form and Integration by Parts Formula 

We start with the rigorous construction of (X!^), which departs from the symmetrizing measure 
Qn on Sjv defined above in (jl.ip . Note that one can identify LP{T,iy,qiy) with LP{T,j\f,qiy), 
p > 1. Throughout the paper denotes both the measure and its Lebesque density, i.e. 
qN{^) = Ja^^^^) measurable A C Sat. 

For all P > 0, N £ N, the measure qi\i satisfies the 'Hamza condition', because it has a strictly 
positive density with locally integrable inverse, c.f. e.g. p.J. This implies that the form £^{f, f) 
with domain / G C°°{T,n) is closable on L'^{T,j\f,qiy). The L^(SAr, qAr)-closure defines a local 
regular Dirichlet form, still denoted by £^ . General Dirichlet form theory asserts the existence 
of a Hunt diffusion (X!^) associated with £^ which can be started in gAr-almost all x G S^r 
and which is understood as a generalized solution of the system ()1.3p . ()1.4p . This identification 
is justified by the fact that any semi-martingale solution solves via Ito's formula the martingale 
problem for the operator {L^ , D{L^)), defined in (II. 2p . 



The following integration by parts formula for gjy can be easily verified (e.g. by approximation 
from corresponding integrals over increasing sub-domains C Sjv). 

Proposition 2.1. Let u G C^C^n) and ^ = be a vector field in C^(En,M.^) 

satisfying {iiv) = on dT,i\f, where v denotes the outward normal field of dTi^. Then, 

N 



{Vu,i)qN[dx) 



div(0 + (lv^-l)E^' 



1=1 



1 



1 



x"- — X 



QNidx). 



Remark 2.2. Let u £ C^{T,j\f) and ^ be a vector field of the the form ^ = w(p with w £ C^(T,j^) 
and ipix) = {(p{x^), . . . , (p{x^)), ip G C°°([0, 1]) and {(p, u) = on dT,]\j, in particular '/3|g[o^i] = 0. 
Then, the integration by parts formula above reads 



{Vu,(,)qN{dx) 



u 



wV^ + {Vw,if) qN{dx), 



where 



(x^ 



/3 



N 



N 



+ Y.^'{x^). 



1=0 i=l 

{/ G C^($]7v) : (Vfju) = on dTj^} as above with v still denoting the outer 
normal field on dT,j\f. Then, for any / G C^{T,j\f) and g G C^^^ we apply the integration by 
parts formula in Proposition 12.11 for ^ = Vg to obtain 



Let C%^^ 



£^if,g) 



fL gq^idx). 



Moreover, it is easy to show that 



In particular, C'i, is contained in the domain of the generator C associated with £^ and 



L^f = Cf for all f £C 



Neu- 
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3 Feller Property 



This section is devoted to the proof of Theorem II. H where we treat the cases (3 < {N + 1) 
and /3 > (A^ + 1) separately. Both cases are not trivial from an analytic perspective due to the 
combination of degenerate coefficients and the Neumann boundary condition. In fact we could 
not find any general result in the PDE literature which contains the current model. Our proof 
for the case /3 < {N + 1) avoids an explicit treatment of the Neumann boundary condition via a 
reflection argument. In the case /? > (A^ + 1) we use a powerful stability property of log-concave 
measures. 



3.1 Case (3 < {N + 1). 

For (3 < {N + 1) the measure qn is no longer log-concave. However, the proof below extends 
also to all cases when f] < 2{N + 1). Let 

Ojv := := Stv n {x G : xjv < 1 - (5}, 

for some positive small 5 < [2{N + 2)]^-*^ and the weight function 

q{x) := qN,s{x) := < ^ ^.^ 3.1) 

l^*^ HiLii^ ^ ^ ) +^ X G Si\Qn, 

where 

Si:={xeR^ : < < < . . . < x^}. 
We want to extend the weight function q to the whole M^. To do this we introduce the mapping 

T:M^^5i (3.2) 

where (1). . . . , (N) denotes the permutation of 1, . . . , such that 

|x«|<...<|x(^)|. 

The extension of q on is now defined via g(x) = q(Tx), x G M^. Again we will also denote 
by q the induced measure on M^. Consider the L^(M^, (7)-closure of 

£'''%fJ)=l {aVf,aVf{x))q{dx), / G C,°°(M^) (3.3) 



still denoted by f"^'" for a measurable field x a{x) G M^^^ on satisfying 

- ■ En <a(xY ■ a(x) <c- En (3.4) 
c 

in the sense of non-negative definite matrices. Let {Yt)t>o = {Y/^'°')t>o be the associated 
symmetric Hunt process on M^, starting from the invariant distribution q. Finally, we denote 
by {Qt)t the transition semigroup of Y. 



5 



3.1.1 Feller Properties of Y 

Let Co(M^) be the space of continuous functions on M.'^ vanishing at infinity. 

Proposition 3.1. For 2{N + 1) > /3 and a matrix a satisfying (|3.4|) . y^'" is a Feller process 
o M^, i.e. 

i) for every t > and every f £ Co(M^) we have Qtf ^ CqCK^), 

ii) for every f G Co(M^), hm^jo Qt/ = / pointwise in M^. 
Moreover, Qtf £ C(M^) for every t > and every f E L'^{R^,q). 

Remark 3.2. It is well known that i) and ii) even imply that limtjo WQtf — f\\oo = for each 
/ G Co(M^). Moreover, the following version of the strong Markov property holds. Let T be a 
stopping time with T < to a.s. for some to > 0. Then, for each / G L^(]R^, q) 

E[f{Yt,)\J^T]=EYAf{Yto-T)], 

with {J-t)t>o denoting the natural filtration of Y . 

Proof of Proposition VJ.1\ ii) follows directly by path-continuity and dominated convergence, i) 
as well as the additional statement follow from the analytic regularity theory of symmetric 
diffusions, see [TH], in particular Theorem 3.5 and Proposition 3.1 and Corollary 4.2, provided 
the following two conditions are fulfilled: 

• The measure q is doubling, i.e. there exists a constant C", such that for all Euclidean balls 
Br C B2R 

q{B2R) < C'q{BR). 

• £ satisfies a uniform local Poincare inequality, i.e. there is a constant C" > such that 

/ \f-{f)Bu\^dq<C'R^ f \Vf\^dq, 

J Br J Br 

for all Euclidean balls Br and / G 'F>{£), where {f)BR denotes the integral /g^ fdq. 

Both conditions are verified once we have proven that the weight function q is contained in the 
Muckenhoupt class A2-, which will be done in Lemma [3^ below. Indeed, the doubling property 
follows immediately from the Muckenhoupt condition (see e.g. [23] or [2T|) and for the proof of 
the Poincare inequality see Theorem 1.5 in l^lOj. □ 

Lemma 3.3. For N such that (3 < 2{N + 1), we have q G A2, i-c. there exists a positive 
constant C = C{N, 6) such that for every Euclidean ball Br 

\Br\ Jbr \Br\ Jbr 
where \Br\ denotes the Lebesgue measure of the ball Br. 

It suffices to prove the Muckenhoupt condition for the weight function q defined by 

q(x) := Ylix' - a;^-i)/3/(^+i)-i, x G Si, (3.5) 
1=1 
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and q{x) := q{Tx) if x G M"^, since there exist positive constants C\ and C2 depending on 5 
and A'' such that 



Ci q{x) < q{x) < C2 q{x), Vx G R^. (3.6) 



Note that (1 — X]\f)N+i ^ is uniformly bounded and bounded away from zero on r^jy. 



Below PR{m) denotes the parallelepiped in with basis point m = (mi, . . . , mjv), which 



is spanned by the vectors Vi = "^jLiej, i = 1,...,N, normalized to the length R, where 



{ej)j=i,...,N is the canonical basis in . Then, PR{m) can also be written as 

PR(m) =m+{xGM^ G [o,^],^^^ + , . . . , G x^-^ + i?] } . 

Lemma 3.4. Let Br he an arbitrary Euclidean hall in M^. Then, there exists a positive 
constant C only depending on N and a parallelepiped Pinirn) C -Si with Z > independent of 
Br such that 

[ q^^dx <C [ q^^dx. 
J Br JPiRim) 

Proof. We denote by (Sn), n = 1, . . . ,2^ ■ N\, the subsets of taking the form 

Sn = {x eM.^ : (si7ri(x), . . . , Siv7riv(x)) G ^i}, 

where Sj G {—1,1}, i = 1,...,N, and 7r(x) = (7ri(x), . . . , 7rjv(x)) is a permutation of the 
components of x. Note that = (Jn and the intersections of the sets Sn have zero Lebesgue 
measure. 

Let now Br{M) be an arbitrary Euclidean ball with radius R centered in M G M-^. We first 
consider the case |M| < 2R. Then, obviously Br{M) C ^4^(0) and B4r{0) n Si C PmiO). 
Hence, 

/ q^^dx < [ q^^dx = V / q^^dx. 



By the definition of q we have 



/ q^^dx= [ q^^dx 

JB4,Rifi)nSn JB4.R{Q)nSl 



for all n G {1,... ,2^ • A^!}. Thus, 

/ q^^dx < 2^ • iV! / q^^dx <C [ q^^dx. 

Suppose now that |M| > 2R. Let no be such that M G Sno- Then, by construction of q we have 

/ q^^dx < / q^^dx 

JBn{M)nSn JBn{M)nSno 

for all n = 1, . . . , 2^ • AT!. Hence, 

/ q^^dx = V / q^^dx < 2^ ■ N\ [ q^^dx. 

Jbu(M) „ JBu(M)nSn JBu(M)nSnQ 
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Set K := T{Br{M) n 5„,)) C Si with T defined as above. Then, it is clear by definition of q 
that 



L 



'Bn(M)nS„Q JK 
and we get 

[ q^^dx <C [ q^^dx. 
JBniM) JK 

Finally, we choose a parallelepiped P2r{t^) such that K C P2R{m) C Si, which completes the 
proof. □ 

Proof of Lemma \3.^A We prove the Muckenhoupt condition for q. Recall that we have assumed 
[5 < 2{N + 1). In the following the symbol C denotes a positive constant depending on N and 
6 with possibly changing its value from one occurence to another. Using Lemma 13.41 we have 

^ ^ qdx [ q-^dx<CR-^^ [ qdx [ q~^dx 



\Br\' 



Br JBr JPrnim) JPmim) 

N 



J PiR(m) j^-, 
N 

" L , ,n< 



i=l 

N 



By the change of variables j/j = — , i = 1, . . . , A^, we obtain 

N 



qdx q-'dx<CR-'''l[ y^'dyi 2/^-+^" 

\i^R\ J Br J Br r=l-^™' 

where we have set fhi := rrii — rui-i with tuq := and := rhi + ^J^-^ ■ for abbreviation. 
Recall that in one dimension the weight function x i— > Ixl*^ on M is contained in A2 rj £ (—1,1) 
(see p. 229 and p. 236 in [21]). Hence, we get for every i £ {1, . . . , N} 

2 



dy^ / y/"^^ ^dy,<c 



IR 



and the result follows. 

□ 



3.1.2 Feller Properties of Y inside a Box 

Let E := := {x G : ||x||oo < 1 - 25} be a box in M centered in the origin and 13{E) be 
the Borel cr-field on E. We denote by 

TE := inf{t >0:YteE^} 

the first exit time of E. This subsection is devoted to the proof of the Feller properties for the 
stopped process Y^^ := Y^'""'^ := Y^'^^, whose transition semigroup is given by 

Qff{x) = E^[f{Yt;,r^)] = E^ [f{Yt)^^r^} + /(y,^)%>,^|] , t > 0, X G ^, 

for every bounded / on E. In order to prove the Feller property we will follow essentially the 
proof of Theorem 13.5 in [7J. It is shown there that the Feller properties are preserved, if the 
domain is regular in the following sense. 
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Proposition 3.5. The domain E is regular, i.e. for every z £ dE we have P^[te = 0] = 1. 

Remark 3.6. z S dE is regular point in the sense of the definition given in Proposition 13.51 if 
and only if for every continuous function / on dE 

lim S-[/(y.Jll|,^<^|]=/(z). 

For the Brownian case we refer to Theorem 2.12 in [T3]. The arguments are robust, but it is 
required that under the probability of the event that the exit time of a ball centered at x 
does not exceed t is arbitrary small uniformly in x for a suitable chosen small t > 0. In our 
situation this property is ensured by [T8|, p. 330. 

The proof of this proposition will be based on the following Wiener test established by Fabes, 
Jerison and Kenig, c.f. Thorem 5.1]. 

Theorem 3.7. Let Bf^^^{0) be a large ball centered in zero such that E C i?^Q/4(0). Then, a 
point z G dE is regular if and only if 

'^J CK£iIJ)f<^' or 

b) /(f" cap(Kp)^(^^ 7 = 

where Kp{z) := {BjiQ{0)\E) n Bp{z) and cap denotes the capacity associated with the Dirichlet 
form S^'"" . 

In the sequel we will use the following notation 

d{x) := |{iG {l,...,iV} ■.x' = x'-^]\, xeM^, 

again with the convention x^ := 0. In order to prove the regularity of E we start with a 
preparing lemma. 

Lemma 3.8. Let y G M.^ be arbitrary. 

i) There exists a positive rg = ro(y) such that for all balls Br{y), < r < tq, there exists 
a parallelepiped Pir{y) contained in Si with d{y) = d(y), a positive constant C and I > 0, 
such that 

q{Br{y)) < CqiPiriy))- 

ii) For all balls Br{y) there exists parallelepiped Pir{y) contained in Si with d{y) = d{y), for 
some I > such that q{Pir{y)) < <i{Br{y))- 

iii) For all y we have Ci r'^(^) < q{Br{y)) for all r and q(Br{y)) < C2 r^^^^ for all r < ro{y) 
for some positive constants Ci and C2 depending on y, where 

h{y) ■.= N-d{y) + j^^d{y). 

Proof. Obviously, due to (|3.6|) it suffices to prove the lemma with q replaced by q defined in 
([33]). 

i) The case y = is clear. For y 7^ we choose tq such that 2ro < \y\ and let no be such that 
y G SuQ. Consider now an arbitrary ball Br{y) with r < rg and let K = T{Sno ^ Bf{y)) be the 
subset of Si constructed in the second part of the proof of Lemma 13.41 Then, possibly after 
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choosing a smaller tq we find a parallelepiped and a positive constant I such that K C Pir{y) C 
S*! with d{y) = d{y) and we obtain i). 

ii) Let K = T{Sna ^ ^riy)) be defined as in i). Then, clearly q{Br{y)) > q{K). Thus, we can 
choose y = Ty and / independent of r such that Pir{y) ^ K and ii) follows. 

iii) We proceed similar as in the proof of Lemma 13.31 For some parallelepiped Pir{y) with 
d{y) = d{y) we use a change of variables to obtain 

JPiAy)) i=i i=i Jyi 

where jji = yi — Vi-i, yo ■= and c, := . Note that d{y) = d{y) is the number of 

components of y which are equal to zero. Using the mean value theorem we obtain that 

for some positive constants Ci and C2 depending on so that iii) follows from i) and ii). □ 

Proof of Proposition |g.5l Let z G dE be fixed and Rq be as in the statement of Theorem 13.71 
Setting h! {z) := 1 — h{z) let us first consider the case h'{z) > —1. Then, using Lemma iii) 
we have that 

/ -TWTT^-^^ s^'(^Us + —-^(Rl-rl)<^, 
Jo QiBs[z)) s Jo 2q{Bro[z)) 

with tq = ro{z) as above in Lemma 13.81 iii). Thus, the criterion a) in Theorem 13.71 applies and 
the regularity of z follows. The case h'{z) < —1 is more difficult. Combining Lemma 3.1 in [lUj 
and Theorem 3.3 in [lOj and using Lemma [3.8l iii) we get the following estimate for the capacity 
of small balls: 

-'"^•<^»Kr^T)"^<r^''S"' '''' 

Recall the definition of the set Kp{z) in Theorem 13. 71 Clearly, for every p sufficiently small there 
exists a ball Bp/2{z) with z depending on p and with d{z) = d{z) such that 5^/2 (^) C Kp[z). 
Let now tq be such that Lemma 13.81 iii) and (j3.7p hold for every ball Bp{z), p < tq. Then, we 
obtain in the case h'{z) < — 1 



f^O ^2 J pro 

/ cap(Kp)— f— - ^ > C / cMBp/2mp^ dp>C 
Jo Q{Bp(z)) p Jo 




/.ro/2 h'{z) 
C{h!{z) + l) . dp 



lo R^^^y^+^ - ph'{y)+l 

Finally, if h'[z) = — 1 we get by an analogous procedure 

Jo ^'-''^''^^WMl^-'' lo P(logi?o-logp)'^ 

/•ro/2 

= C / (-log(logi?o - logp)) dp = 00. 
Jo 

Hence, applying the criterion b) of Theorem 13.71 completes the proof. □ 
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Proposition 3.9. is a Feller process, i.e. 

i) for every t > and every f G C{E) we have Qff G C{E), 

a) for every f G C{E), lim^jo Qf f = f pointwise in E. 

The statement is classical and can be found e.g. in Theorem 13.5 in [7]. For readability we 
repeat the argument here. We shall need the following lemma. 

Lemma 3.10. For any compact set K <Z E we have 

limsup P^'Ite < i] = 0. 

^iO x(^K 

Proof. We need to show that for any (5 > there exists a to > such that 

inf P^[Ti5 > to] > 1 -5- (3.8) 



Consider a bounded function / G Co(M ) such that 0</<l,/ = loni^ and / = on the 
complement of E. Let now to be such that sup^<^Q WQtf — f\\oo < (cf. Remark l3.2p . Then, 

E^ifiYtJ] < P^[te > to] + E^ifiYtJ 

For X €z K the left hand side is equal to 

QtJix) = 1 + QtJix) - fix) > 1 - sup WQtf - /lloo > 1 - t 

t<to 

On the other hand, using the strong Markov property (cf. again Remark 13. 2p and the fact that 
f{Yl,) = we have 



E^ifiYt,) %^<,„|] = E^ [E^ifiYt,) %,<i„|] = E^ [Qt,^,J{Yr,) 1 



E^ [{Qt,.rJ{Yr,) - f{Yr,)) t{r,<to}] < 



{TE<to}\ 



sup IIQt/- /lloo ]l{rs<to} 

t<to 



6 

and (BIHD follows. □ 



Proof of Proposition [XPl Let t > and / G C{E). Then, by the semigroup property of (Qf) 
we have for < s < t 

Qff{x)=E- [i,s{Ys^rE)], 

where 

M^) = Qlsfix) = E- [f{Y(^t-s)Arj] , XeE. 

Then, ipg can be extended to a function in L^(M^,g) and by Proposition 13.11 we have Qs^^s £ 
a(M^). Since 



\Qffix) - QsMx)\ = lE^^lMYsAr^) - Mys)]\ < 2 IIV.II P^'Ite <s]<2 II/II P^'It^ < s] 

and since the right hand side converges to zero uniformly in x on every compact subset of E 
by Lemma [3.101 we conclude that Qff G Ci,{E), i.e. Qff is continuous in the interior of E. 
In order to show i) it suffices to verify continuity at the boundary. Since E is regular, we have 
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obviously Qf f = / on dE. By Lemma 13.1 in ff] we have upper semicontinuity of the mapping 
X 1-^ P^[t < te]- Hence, we obtain for z £ dE, 

nmsupP'^[t < te] <P''[t< te] = 0, 

where we have used the regularity of E in Proposition 13.51 Thus, for every x G E 

\Qff{x) - f{z)\ < \Qff{x) - i?-[/(y,J%^<„o|]| + \E-'[f{Y^,)l{,^^^y] - f{z)\ 

< \E^ [/(yt)]i|,<,,| + fiYr,)i{t>r,y - /(y.j%,<„o|] I 

+ \E-[f{Yr,)l{,^^^}]-f{z)\ 

<||/IUP"[t < te] + \E^ [/(y.jii|.,<o,| (i|i>.,| - 1)] I 

<2||/|UP"[t < te] + \E-[fiYr,)l{r,^oo}] - fiz)\ , 

where we have used the fact that the event {te < oo} is included in {t > te}- Since z is a 
regular point, the second term tends to zero as x ^ z, x £ E, cf. Remark 13.61 Hence, 

hm Qf/(x) = /(z) 

and property i) is proven. 

To prove ii), we extend / to a function in Co(M^), i.e. we may deduce from Proposition 13.11 
that limtjo Qtf{x) = f{x) for every x £ E. Furthermore, we have for every x £ E 

\Qff{x)-Qtf{x)\<P^[t<TE]\\f\U 

and P^[te = 0] = 0, since E is open and Y has continuous paths. Hence, 

limQf/(x) = limQJ(x) = /(x). 

This gives pointwise convergence on E. Since Qf f = / on dE by regularity, the convergence 
on dE is trivial. □ 



3.1.3 Feller Properties of 

In this section we will finally prove the Feller property for stated in Theorem ll.il To this 
end we construct a Feller process X taking values in Eat and in a second step we will identify 
this process with X'^ . 

In analogy to the definition of fijv above we set 

Qi:=nf ■.= {x e"EN ■■ x'+^ -x' >6}, i = 0,...,iV. (3.9) 

and moreover 

A' := dnf\dJ:N = {xe^N ■■ -x' = 26}. 

Notice that we can choose 6 so small that T,^ = IJiLi ■ Furthermore, we define the mappings 
H„i = 0,...,N, by 

Hi{x) := {x\ ...,x\l-{x^ - x'), 1 - (x^-i - x'), . . . , 1 - {x'+^ - x')) , x G S^. 

Notice that for every i, Hi maps Qn on Qj and vice versa. In particular. Hi o Hi = id and 
Hi\i is the identity on JIat. Let T : E ^ C f^Tv be defined as in (13. 2p . Let y* denote the 
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R^-valued Feller process induced from the form ()3.3p with a = ai, where the matrix valued 
function a is defined g-almost everywhere by 

J DH^ for X e 5i 

MX) :- I . for x G 5„, n e {2, • • • 2^ • iV!}, ^^"^^^ 

such that condition (|3.4|) is clearly satisfied (note that DHi = DH^^). In particular, Oj is 
constant on every Sn- Moreover, setting pi := Hi oT, clearly Oj = (Dp^^y. 

Remark 3.11. Analogously to Proposition 12.11 one can establish the following integration by 
parts formula associated to E^'""^. Let u E C^(M^) and ^ a continuous vector field on such 
that supp^ ^ {x G : ||2;||oo < 1 — (5}, ^ is continuously differentiable in the interior of each 
Sn and for every n ^ satisfies the boundary condition (a* • Oj • Ct'n) = on dSn, I'n denoting 
the outward normal field of dSn- Then, 

/ {aiVu,ai^)q{dx) = - u [div(a- • • ^) + (a- • Oj • Vlogg)] 

JRN J^N 

Thus, every smooth function g such that ^ = Vg satisfies the above conditions is contained in 
the domain of the generator L* associated to S'^^^-i and on \ dSn we have 

Vg = div(a* • at ■ Vg) + (a* • ai ■ Vg, V log q). 

Next we define the r2?''-valued process X* by 

Xl = pi{Y:''^) = HioT{Y:'''), t>0. 

The semigroup of X' wiU be denoted by (P/)t>o, i-e. Pif{x) = 
Lemma 3.12. For every i, is Markovian. 

Proof. Since Hi is an injective mapping for every i, it suffices to show that the process T{Y.^'^) 
is Markovian. Moreover, since T{Y^'^) = T{Y^\^^) = {T{Y.'^))t/\TE it is enough to prove the 
Markov property for the process r(y*), which is implied e.g. by the condition that for any 
Borel set A C Sat 

px^Yi G T-\A)] = py[Yi £ T-^{A)] whenever T{x) = T{y). (3.11) 

Now the choice of q and the metric Oj imply for any Borel set A Q Y^n condition (|3.1ip is 
satisfied. To see this, let {ak \ k = 1, - ■ ■ ,N + N{N — 1) /2} be the collection of line-reflections 
in M.^ with respect to either one of the coordinate axes {Ae^, A G M} or a diagonal {A(ej ■+ e^)}, 
then for x, y G with T{x) = T{y) there exists a finite sequence such that r(x) := 

Cfci ° o'/c2 ■ ■ ■ ° <^ki{x) = y- Now each of the reflections ai preserves the Dirichlet form (j3.3p 
when a is chosen as in (I3.10p . such that the processes r(Y!*'^) and Y!*'^ are equal in distribution. 
Moreover, t{T~^[A)) = T~^(^), from which (j3.1ip is obtained. □ 

Lemma 3.13. For each f G C^g„ the process t f{Xl) — Jq f{Xl) ds is a martingale w.r.t. 
to the filtration generated by X\ 

Proof. Similar to Proposition 12.11 one checks for / G C^^^ H Cc{^i) that the function /j on M^, 
which is defined hy fi = f o Hi o T = f o pi on the set {x G : ||x||oo < 1 — 5} and fi = 
on the complement this set, belongs to the domain of the generator U of the Dirichlet form 
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(|3.3p with a = ai as in p.lOp (cf. Remark [3. lip . Hence the process fi{Y.^) — ^QUfi{Yg)ds is 
a martingale w.r.t. to the filtration generated by and thus also fi{Y.^'^) — U fi{Yp^)ds 
due to the optional sampling theorem. Obviously in the last statement the function / can 
be modified outside of Qf^, i.e. it holds also for / G Cnbu- Moreover, fi{Y.^'^) = f{X!) and 
Pfi = (L^f) o Pi on E. Thus, /(X) - Jq f{Xi)ds is a martingale w.r.t. the filtration 
generated by Y^ and adapted to the filtration generated by X* = HioT{Y.^'^) which establishes 
the claim. □ 

Proposition 3.14. For every i, X* is a Feller process, more precisely 

i) for every t > and every f £ C{0,f^) we have P^f £ C{0,f^), 

a) for every f G C{il.f^), lim^j^o^t/ = / pointwise in 0?^. 

Proof. Since obviously for every continuous / on Qf^ 

PI fix) = E^[f{xi)] = E-[f{H, o r(y/'^))] = Qf{f o o r)(x), t > o, 

the result follows from Proposition 13.91 and the continuity of the mappings Hi and T. □ 

Next we define the process X with state space Sjv as follows: Let be the initial distri- 
bution of X, i.e. Xo ~ QN- Choose ii G {0, . . . X} such that Xq G and dist(Xo,^*i) = 
maxj dist(Xo, A*). We set Xf = Xl^ for < t < Ti, where Ti denotes the first hitting time of 
i.e. on [0, Ti] the process behaves according to P*^. Choose now i2 £ {0, . . . X} such that 
Xt^ G 0?^ and dist(X7^j, A*^) = maxj dist(XT^ , ^*). At time Ti the process starts afresh from 
Xt^ according to (P/^) up to the first time T2 after Ti, when it hits A^'^ . This procedure is 
repeated forever. 

Proposition 3.15. X is a Feller process. 

Proof. Let {Pt) denote the semigroup associated to X. For / G C(Siv), we need to show that 
Ptf G C(SAr) for every t > 0. Let us first show that Pt„/ £ C(S7v) for every n using an 
induction argument. For an arbitrary x G Sat, choose ii as above depending on x such that 
PT^f{x) = Py^/(x). Since P^^f G (7(0?^^), we conclude that Pti/ is continuous in x for every 
X. For arbitrary n and x G Sa? we have by the strong Markov property 

PTr.^J{x) = E^[f{Xr^^,)]=E^ 

and since ^r"+i-T„/ '^^'^ extended to a continuous function on Sat, we get Pr^+if £ C(Sjv) 
by the induction assumption. 

Similarly, one can show that for every n the mapping x 1-^ E^[f{Xt)l{te(T„,Tn+i]}] is continuous. 
Finally, for every x G S^r 

n-l 

Ptf{x)-J2E- 

k=0 

and since r„ ^ cxd P^'-a.s. locally uniformly in x as n tends to infinity, the claim follows. □ 

The proof of Theorem ll.ll is complete, once we have shown that the processes X and the original 
process X^ have the same law. 



E^ 



P^r 



/(X)]l{te(Tfe,Tfe+i]} 



< 
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Lemma 3.16. For x G Sjv the process (X^) obtained from conditioning X to start in x, solves 
the martingale problem for the operator (L^,C^g^) as in f|1.2p and starting distribution 6x- 

Proof. Let {Tf^} be the sequence of strictly increasing stopping times introduced in the con- 
struction of the process X. Then for s < t 

^ {Tk+iVs)M 

f{Xt) - f{Xs) - ^ L^f{X^)da = ^ /(l(T,+,vs)At) - /(^(T,vs)At) - / L^f{X„)dcT. 



{TkVs)At 



Hence 



E{f{Xt) - I L''f{X^)da \j^s) = f{Xs) - f L''f{X„)da 
Jo Jo 

/ L''f{X^)da 



(TfeVs)At 

Using the strong Markov property of the Feher process X. and its pathwise decomposition into 
pieces of {X*}-trajectories one obtains 

{Tfe+iVs)At 



K(/(^{T.+iVs)At) - /(%,vs)At) - I L''f{X^)da \Ts) 

(TfeVs)Ai 

=^[%(.,v.)..(/(^:A(t-.))-/w)- r^*''^^^/w)^^) (3.12) 

J 



where, by construction of X. , r ist the hitting time of the the set j4j^, for which dist(^jj. , Xj^^) = 
maxj dist(j4j,Xrj,). Lemma 13.131 imphes that the inner expectation in (j3.12p is zero. □ 

The last ingredient for our proof of Theorem 1 1.1 1 is the identification of the processes X. with X.. 
Since both are Markovian and solve the martingale problem for the operator (L^, C^g^) it suf- 
fices to show that the martingale problem admits at most one Markovian solution. Clearly, any 
such solution induces a symmetric sub-Markovian semigroup on L^(S7v,Q'7v) whose generator 
extends (L^,C^g^). Hence it is enough to establish the following so-called Markov uniqueness 
property of (L^,C^g„), cf. [S*, Definition 1.2], stated in Proposition 11.21 

Proof of Proposition [TT^ Let //^'^(Sjv, (/Af) (resp. 'ffQ'^(SAr, qtv)' in the notation of [8]) denote 
the closure of C%^^ w.r.t. to the norm \\f\\, = (||/|li2(s^,,^) + ||V/||i2(s^,g^))^/2 

W^''^ be the Hilbert space of L^($]jv, gAr)-functions / whose distributional derivative Df is in 
L^iT,N,qN), equipped with the norm \\f\\^ = (||/||i2(2^^^^) + ll-D/||^2(s^,g^))^/^- Clearly, the 
quadratic form Q{f,f) = Z)/)^2(5^^ g^), 'D{Q) = W^''^{T,iy,qj\f), is a Dirichlet form on 
L^(SAr,qjv)- Hence we may use the basic criterion for Markov uniqueness [HI Corollary 3.2], 
according to which (L, C^^^) is Markov unique if H^'^ = W^''^. 

To prove the latter it obviously suffices to prove that H^'^ is dense in T^^'^. Again we proceed by 
localization as follows. For fixed 6 > OletQf = J7j denote the subsets defined in ()3.9p . then {f^f*^} 
constitutes a relatively open covering of S^v for 5 small enough. Let {r]i}i=Q^... ^tv and {xi}i=o,- ,N 
be two smooth partitions of unity on S^v such that r/j = 1 on Qf^ and supp(7/j) C ftf^ and Xi = ^ 
on rj^^ and supp(xi) C respectively. Hence writing / G W^''^{Y,i\f,qi\f) as / = /i with 
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fi = Vi' f ^ qj\f) it suffices to prove tliat fi G H^''^(T,]\f, q^). We sliow first that fi can 

be approximated w.r.t. ||.||-|^ by functions whicli are smooth up to boundary and secondly that 
such functions can be approximated in ||.||^ by smooth Neumann functions. We give details for 
the case i = N only, the other cases can be treated almost the same way by using the maps Hi. 

Let gN = f ■ Xn G W^^'^l^AT, gAr), then the restriction of gN to fi^ belongs to the space 
qat) = W^''^{Qf^,q), where q denotes the modification of q^ according to (j3.ip . 
Due to Lemma 13.31 the extension q{x) = g(T(x)), x G M^, lies in the Muckenhoupt class 
A2. Further note that W^^^{n%,q) C W^^"^ {n% , dx) = H^^^{n%,dx) if p < N + 1, and 
c W^^^{n%,dx) = H^^^{Q%,dx) by the Holder inequality if + 1 < /3 < 2(iV + l), 
such that gN has well defined boundary values in L^{dQ,N,dx). Hence we may conclude that 
the extension gN{x) = gN{T{X)), x G defines a weakly differentiable function on with 
||ffAf|liyi.2(]R^ q) ' IIS'A^II wi.2{Qjv !?)■ Theorem 2.5] the mollification with the stan- 

dard mollifier yields an approximating sequence {ui}i of smooth functions ui G C^(R^) of g^ 
in the weighted Sobolev spaces i7^'^(M^, q). Extending each ui-rj^ by zero in Sjv\$7^ we obtain 
a sequence of C°°(S7v)-functions {ui • r/jy}/ which converges to gN • Vn = In in H^''^{'^n,Qn)- 
This finishes the first step. In the second step we thus may assume w.l.o.g. that g^ is smooth 
up to the boundary of T,n- In particular, g^ is globally Lipschitz. Since qN is integrable on 
Sjv we may modify /tv close to the boundary to obtain a Lipschitz function /n which satisfies 
the Neumann condition and which is close to /at in ||.||-^. (Take, e.g. /Ar(x) = /(7r(x)), where 
7r(x) is the projection of x into the set = {x G T,n, |dist(x, OTin) > e} for small e > 0.) We 
may now proceed as in step one to obtain an approximation of /at by smooth functions w.r.t. 
\\.\\i, where we note that neither extension by reflection through the map T nor the standard 
mollification in [T3] of the extended /at destroys the Neuman boundary condition. □ 

Hence we arrive at the following statement which implies the first statement of Theorem 11.11 for 
the cases when /? < 2{N + 1). 

Corollary 3.17. For quasi-every x G Sat, the the processes and are equal in law. In 
particular, X is a Feller process on Tin- 

3.2 The case /9 > (iV + 1). 

The estimate (|1.5|) looks (also for the case (3 < + 1) like a straightforward application of 
the Bakry-Emery r2-calculus. However, the complete Bakry-Emery criterion requires the 
condition on an algebra of functions which is dense in the domain of the generator C of £ w.r.t. 
the graph norm. The verification of the latter typically leads back to elliptic regularity theory 
for C which we want to avoid. Instead our argument below is based on a recent powerful result 
by Ambrosio, Savare and Zambotti [2] on the stability of reversible processes with log-concave 
invariant measures. 

Proof of (jl.Sp . By abuse of notation let qN be the extension of the measure qN to given by 

qN{A) = qN{A n S^) = ^ / e'^^-^^^^^'^dx 

JAn{V<oo} 

for any Borel set yl C M^, where ^ : ^ M U {oo}, 

V{x) = 



Eill^ ln(xi - Xi-i) X G 



oo xeR^\T.N 
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is a lower semicontinuous and convex function on M^. In fact, by elementary calculations, for 
X G Sat with pi = pi{x) = 1/ (xj - Xi^if, i = 1, . . . , + 1, 



N-l 



i=l 



where /cat > denotes the smallest eigenvalue of the strictly elliptic matrix A = (2(5ij — G 
I^NxN_ Hence the measure qn is log-concave in the sense of [2]. (Note that according to ibid. 
Theorem 1.2.b. the Feller property of X. on S^r is automatically implied.) Let x £ Ti^ denote 
the barycenter of the simplex S^r. For e > let 



V,{x) 



V{{l-e)x + ex)) X G Sat 
oo X G \ S 



N 



and define the measure on by 



Since K G C°°(SAr) and the boundary of Sjv is piecewise smooth this process can also be 
constructed by the corresponding Skorokhod SDE. Moreover, the classical coupling be reflection 
method can be applied, c.f. [25J. Taking expectations in [23 eq. (2.5)] (note that / = in 
our case) and using the strict convexity of together with Gronwall's lemma this yields the 
estimate 

Ec(Z:'\Zf)<e-(^-^)(^-^)^^-V(x,y), 

where C denotes the law of the coupling process {X^'^ , Xl''^)t>o, starting in (x, y). In particular 
the following estimate in the L^-Wasserstein distance di for the heat kernel of (X^) is obtained 

di(P/(x,.),Pi^(y,-)) < e-(i^-i)(^-^)''=-*d(x,y). 

Since — > qn weakly, we may now invoke Theorem 6.1. of [2] in order to pass to the limit for 
e ^ in the left hand side above, using also the continuity of the L^- w.r.t the L^-Wassserstein 
metric. Hence we arrive at 

di(Pi(x,-),Pt(2/,-))<e-(i^-^)(^-^)'^'^*d(x,y), 
for all e > small enough, thus also for e = 0. Via Kantorovich duality this implies 

\Ptf{x)-Ptf{y)\ < e-(i^-^)'=^*Lip(/)d(x,y), 
for or all / G Lip(i;7v) and x,y £ Eatj which is the claim. □ 



4 Semi-Martingale Properties 

In this final section we prove the semi-martingale properties of X^ stated in Theorem 11.41 To 
that aim we establish the semi-martingale properties for the symmetric process X^ started in 
equilibrium, which imply the semi-martingale properties to hold for quasi-every starting point 
X G Sat and by the Feller properties proven in the last section for every starting point x G Sa^- 
In order to establish the semi-martingale properties of the stationary process, we shall use the 
following criterion established by Fukushima in [13]. For every open set G C Sat we set 

Cg := {u G D{S^) n C(S^) : supp(^/) C G}. 
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Theorem 4.1. For u E D{£^) the additive functional u{X^) — u{Xq) is a semi-martingale 
if and only if one of the following (equivalent) conditions holds: 

i) For any relatively compact open set G cTi^, there is a positive constant Cq such that 

\£''{u,v)\ <Cg hWoo. yveCc. (4.1) 



a) There exists a signed Radon measure v on S^v charging no set of zero capacity such that 



(4.2) 



Proof. See Theorem 6.3 in [13J. 



□ 



Theorem 4.2. Let be a symmetric diffusion onTi^ associated with the Dirichlet form 
£^ , then X^ is an M.^ -valued semi-martingale if and only if f3/{N + 1) > 1. 

Proof. Since the semi-martingale property for M^-valued diffusions is defined componentwise, 
we shah apply Fukushima's criterion for u{x) = x^ , i = 1, . . . , N . 

Let us first consider the case where /3' := f3/(N + 1) > 1. Then, for a relatively compact open 
set G C Sat and v G Cq, 



_d_ 



v{x) q^idx) 



^ / dx' dx'--- / dx'+i / / dx"" r\ { 



i=o 



dx^ ^ ^ ^ ' 

Since /3' > 1, we obtain by integration by parts 

jv{x) {x' - x'^y-'^ (x^+1 - xY~^ dx' 



i-i dx^ 



r.i+1 



{13' - 1) / v{x) {x' - x'^-^ (x*+i - xY^^ - (.x' - x'~y-^ {x'+^ - xY'^ 



dx' 



so that 



d 



^t;(x) {x' - x'-y {x 



-xY-'dx' 



<w'-mv\ 



{x'-x'-Y~^dx'+ I {x 



r-i+^ - xY-^ dx' 



<2(/3'-1)||7;||oo / r^'-^dr < C \\v\\oo, 
Jo 

and we obtain that condition ()4.ip holds. For /?' = 1 the measure qn coincides with the 
normalized Lebesgue measure on E^r, condition (j4.ip follows directly. 
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Let now /?' < 1 and let us assume that u{X^) —u{Xq) is a semi-martingale. Then, there exists 
a signed Radon measure v on Sat satisfying (j4.2p . Let = z^i — 1/2 be the Jordan decomposition 
of V, i.e. vi and V2 are positive Radon measures. By the abobe calculations we have for each 
relatively compact open set G C E^r and for all v G Cq 



1 r ^ 

£^{u,v)=-—{p'-i) / v{x) n 



dx. 



j=0 

Hence, we obtain for the Jordan decomposition = 1^1 — 1^2 that 

1 r ^ 

„^{G) = —(1 - /?') / (x^+i - xY~^ TT(^^'^^ - a^^y^^f^^; 
Jg fJ-Q 

1 /• ^ 

:.2(G) = —(1 - (3') / (x^ - x*-i)^'-2 rr (x^+i - x^y-^dx. 

Set := {x G ^Sat : x^ = x-'^-'^}, j = 0,. . . ,N, and let for some xq G 9S* and r > 0, 

A := xo + [—r, r]^ n S^v be such that dist(^, > for all j 7^ i. Furthermore, let (j4„)„ be 

a sequence of compact subsets of A such that A^ ] A and dist(A„, dT,]^) > for every n. 
By the inner regularity of the Radon measures vi and 1^2 we have vi{A) = lim„i/i(yl„) and 
V2{A) = lim„z^2(^n)- Since j3' — 2 < —1, we get by the choice of A that fiiA) = 00, while 
i'2{A) < 00, which contradicts the local finiteness of v and z^i, respectively. □ 
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